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Abstract 

> . 

£^ , Nucleation rates for tunneling processes in Minkowski and de Sit- 

ter space are investigated, taking into account one loop prefactors. In 
qq ■ particular, we consider the creation of membranes by an antisymmet- 

ric tensor field, analogous to Schwinger pair production. This can be 
viewed as a model for the decay of false (or true) vacuum at zero tem- 
perature, in the thin wall limit. Also considered is the spontaneous nu- 
cleation of strings, domain walls and monopoles during inflation. The 
instantons for these processes are spherical world-sheets or world-lines 
embedded in flat or de Sitter backgrounds. We find the contribution of 
. such instantons to the semiclassical partition function, including the 

one loop corrections due to small fluctuations around the spherical 
worldsheet. We suggest a prescription for obtaining, from the parti- 
tion function, the distribution of objects nucleated during inflation. 
This can be seen as an extension of the usual formula, valid in flat 
space, according to which the nucleation rate is twice the imaginary 
part of the free energy. In addition, we use the method of Bogolubov 
transformations to compute the rate of pair production by an electric 
field in 1+1 dimensional de Sitter space, and compare the results to 
those obtained using the instanton method. Both results agree where 
they are expected to, not only in the exponential dependence but also 
in the prefactor, confirming the validity of instanton techniques in de 
Sitter space. Throughout the paper, both the gravitational field and 
the antisymmetric tensor field are assumed external. 
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1 Introduction 



A wide class of non-perturbative phenomena in field theory can be under- 
stood in terms of quantum tunneling. A well known example is the decay 
of false vacuum: the materialization of bubbles of true vacuum in first order 
phase transitions |], |!|. Lower dimensional versions of this process have been 
used to model the decay of metastable topological defects, such as domain 
walls and strings 0. 

A closely related phenomenon is the neutralization of the cosmological 
constant through membrane creation. In a spacetime of dimension d = N+l, 
an antisymmetric tensor field A of rank N induces a cosmological constant. 
This is because the corresponding field strength F = dA has only one inde- 
pendent component, which has to be constant in the absence of sources. Just 
as an electric field decays through Schwinger pair creation, this cosmological 
constant decays through membrane creation if A is coupled to a membrane, 
a process first described by Brown and Teitelboim M . 

Such tunneling processes can happen in flat as well as in curved spacetime. 
In addition, in curved spacetime, new effects can arise. It has been shown that 
topological defects such as circular loops of string, spherical domain walls and 
monopole-antimonopole pairs can spontaneously nucleate during inflation in 
the early universe 0. These nucleations are somewhat analogous to particle 
production by an external gravitational field. Another consequence of space- 
time curvature is the possibility of true vacuum decay ||, through nucleation 
of false vacuum bubbles. 

Nucleation processes can be described using the instanton methods M. 
The instantons are classical solutions of the Euclidean equations of motion, 
with appropriate boundary conditions. They are saddle points of the Eu- 
clidean path integral, and as such they provide the basis for a semiclassical 
evaluation of the partition function. The contribution of one instanton to 
the path integral has the form 

Ae~ Se , (1) 

where Se is the Euclidean action of the instanton, and the prefactor A arises 
from Gaussian integration over small fluctuations around the instanton. The 
main part of this paper will be devoted to the calculation of the prefactors 
A for the class of processes mentioned above. 
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The instanton methods can be applied in flat and in curved backgrounds. 
One limitation of the formalism is, however, that the spacetime under consid- 
eration has to have real Euclidean sections. Here, we shall concentrate on de 
Sitter and Minkowski space. One reason for studying de Sitter space is that 
it describes the geometry of spacetime during inflation. A de Sitter space of 
dimension d can be defined as a hyperboloid embedded in a Minkowski space 
of dimension d + 1 , 

VAB X A X B = H- 2 . (2) 

Here X A are the coordinates in the embedding Minkowski space (A = 0, d), 
H is the expansion rate during inflation and tjab is the Minkowski metric. 
We use the metric convention (—,+,...,+). The Euclidean section of de Sit- 
ter space is obtained by analytically continuing the temporal coordinate X° 
to imaginary values 

X° = -iX° E 

(with X% a real number). With this rotation the hyperboloid (Q) becomes a 
d-sphere of radius if -1 embedded in flat Euclidean space (see Fig. 1). 

In flat space, at zero or at finite temperature, the nucleation rates can be 
related to the imaginary part of the free energy |7], and they are essen- 
tially given by an expression of the form ([]]). In curved space, it is believed 
that the nucleation rates also have the exponential dependence ([!]), although 
the theory has not been developed to the same level of rigor than in flat 
space. When the size of the instantons is very small compared with 
one expects that the usual flat space formulas should apply. However, for the 
spontaneous nucleation of topological defects and for the nucleation of false 
vacuum bubbles, the size of the instanton is comparable to the horizon scale 
H~ x . In such cases it is not clear how one should compute the nucleation 
rate, and one may even question whether such tunnelings can occur. 

In this paper we take the heuristic point of view that these nucleations can 
indeed occur. We suggest a prescription for computing, from the semiclas- 
sical partition function, the equilibrium distribution of created membranes, 
bubbles and topological defects during inflation. This can be seen as a gen- 
eralization of the formulas that one uses in flat space. 

In 1+1 dimensions the process of membrane (or bubble) creation reduces 
to that of particle creation in an external field. In that case, the predictions 
of the instanton method can be compared with the results that one obtains 
by using the better understood method of Bogolubov transformations. As 
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we shall see, the results of both methods agree (in the limit where they are 
valid), even in the case when the size of the instanton is comparable to the 
horizon scale. Then, at least in 1+1 dimensions, the instanton prescription 
seems to be valid, and there is no reason to believe that it will not be valid 
in higher dimensions. 

Finally, we should mention that throughout the paper, both the gravita- 
tional and the antisymmetric tensor fields are asumed to be external. As we 
shall see, the backreaction of the membrane on the antisymmetric field can 
be neglected when the charge e of the membrane is very small compared to 
the field strength. In the context of false vacuum decay, this means that the 
difference in energy density between the vacua in the two sides of the mem- 
brane will be much smaller than the overall cosmological constant. Also, in 
order that gravitational field of the membrane be negligible, the mass scale 
of the membrane should be sufficiently small (see e.g || for a comparison of 
instantons with and without self gravity). Gravitational backreaction effects 
are interesting in their own right, leading to qualitatively different behaviour 
at large mass scales, but introducing new complications and problems 0, §]. 
These will be left as subject for future research. 

The plan of the paper is the following. In section 2 we study the instantons 
for membrane production. These are the same as the ones studied in Ref. 
[B], in the limit of negligible self gravity and small charge. These instantons 
are essentially spherical Euclidean worldsheets of dimension N and radius Ro 
(representing the membrane) embedded in Euclidean de Sitter space, which 
is itself a sphere of radius if -1 and dimensionality d — N+l. The radius Ro 
is determined by the strength of the antisymmetric field, the charge of the 
membrane, its surface tension and the expansion rate H. Special attention is 
paid to the analytic continuation of the instantons to Lorentzian signature, 
and the effect of de Sitter transformations on the resulting solutions. These 
Lorentzian solutions describe the motion of the membranes after nucleation. 

As mentioned above, to calculate the prefactor A we need to integrate 
over small fluctuations around the instanton. In Section 3 the theory of such 
fluctuations is reviewed, with an emphasis in the so-called zero modes. The 
zero modes are perturbations which do not change the shape of the instanton, 
but correspond to infinitessimal translations of the solution as a whole. We 
make use of the covariant formalism developed in Refs. [JTD, IT], |12| , according 
to which the worldsheet fluctuations are represented by a scalar field 'living' 
in the unperturbed worldsheet. This scalar field has the meaning of a normal 
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displacement of the worldsheet. 

In Section 4, we briefly review the instantons for the spontaneous nucle- 
ation of topological defects during inflation. These can be seen as a limiting 
case of the ones for membrane creation, when the external antisymmetric 
field is switched off. However, in this case the co-dimension of the defect's 
worldsheet can be larger than one and additional perturbations have to be 
considered. For later use in the paper, we also discuss the instantons for a 
massive particle at finite temperature and for pair production in 3+1 dimen- 
sions. We study fluctuations around these instantons and give the normal- 
ization of the zero modes. 

Section 5 is devoted to the semiclassical evaluation of the partition func- 
tion. It is shown that the evaluation of the prefactor is formally equivalent 
to the evaluation of the effective action for a free scalar field in curved space- 
time (this curved space-time is the worldsheet of the instanton; a sphere, 
in our case). To illustrate the method, the instanton formalism is used to 
compute the partition function for a gas of massive particles at finite tem- 
perature, in the semiclassical limit M. » T, where M. is the particle's mass 
and T is the temperature. 

In section 6 we compute nucleation rates in flat space, recovering known 
results for pair creation in 1+1 and 3+1 dimensions, and for bubble forma- 
tion in 2+1 and 3+1 dimensions. The question of renormalization is briefly 
discussed in analogy with the renormalization of the effective action for a 
scalar field in curved space. 

In section 7 we discuss the nucleation rates in curved spacetime. We find 
the size distribution of created membranes, bubbles and defects during infla- 
tion. We also give the momentum distribution for the case of pair creation. 

In section 8 we study the quantization of a charged scalar field interacting 
with an external electric field in 1+1 dimensional de Sitter space. The spec- 
trum of particles created by the electric and gravitational field is computed 
using the method of Bogolubov transformations. 

Some conclusions are summarized in Section 9. The computation of func- 
tional determinants on the N-sphere, necessary for the evaluation of the in- 
stanton prefactors, is done in the Appendix. 
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2 Production of membranes by an antisym- 
metric tensor field. 



In this section we describe the instantons for the creation of membranes by 
an antisymmetric tensor field and their analytic continuation to Lorentzian 



signature. The instantons discussed in subsection 2.3 also represent the for- 



mation of true (or false) vacuum bubbles in the thin wall limit. 

2.1 Particle coupled to an electric field in 1+1 dimen- 
sions 

For the main part of the paper, the antisymmetric tensor field will be assumed 
external. However, it is instructive to start our discussion with the 1+1 
dimensional case and treating the electric field as dynamical. The action for 
a spinless particle of mass M. an charge e interacting with a Maxwell field 
in 1+1 dimensions is given by M 

S = -M J ds + ej A^dx" -^J d 2 x^F^ u + j d 2 xd pL {^F^ A u ). 

(3) 

Here x^(s) is the particle's trajectory and ds is the proper time interval. The 
particle's world-line is indicated by E and g denotes the determinant of the 
spacetime metric g^ v . The last term is a boundary term, included to yield 
a well defined equation for A^. (This is necessary since we shall consider 
situations in which the field strength F^ u = d^A u — d u A^ does not vanish at 
infinity Q.) 

Using the notation of forms, A = A^dx 11 , the field strength is F = dA. 
In 1+1 dimensions F has only one independent component, the electric field 
E, 

F = -El (4) 

where 

e = yjjg'le^dx' 1 A dx u . (5) 

Eq.(]|D is just saying that in 1+1 dimensions any 2-form is proportional to 
the "volume" form e. (Here is the antisymmetric symbol, with eoi = 1.) 
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Variation of (|3]) with respect to yields the Maxwell equation 

d,(^gFn = -eJ ds5 2 (x - x{s))^±. (6) 

This equaton implies that in the absence of sources the electric field is con- 
stant 

E = const. 

whereas the effect of a charged particle is to produce a discontinuity 

AE = e 

in the electric field as we cross from one side of the worldline to the other. 
Therefore, the only equation for A^ is Gauss' law, which is an equation of 
constraint |13| . In this sense A^ does not actually have any field degrees of 



freedom in 1+1 dimensions. 

As is well known, a constant electric field will decay by producing pairs of 
charged particles, which will then accelerate away from each other, screening 
the electric field away. The rate of pair production can be computed using the 
instanton methods. For this, it is necessary to solve the Euclidean equations 
of motion. The Euclidean action Se can be found by complexifying the 
temporal coordinate x° — > —ix° E , ds — > —idsE and leaving the field strength 
F^ v — > Fg V unchanged Q. From F 01 = d°A 1 — d 1 A° one has to complexify 
the vector potential A — > A° E , A 1 — > iA\. With such rotations the Euclidean 
action is found to be {iS — > —Se), 

S e = +M J ds + ej A^dx» - l -J d 2 ^F, u F^ + J dx 2 d,(^F^ A v ). 

(7) 

Here we have dropped the Euclidean subscripts. Expanding the derivative 
in the last term and using (pi) one has 



S E = Mj ds + ^J d 2 x^F, u F^, 



which is positive definite. 

Unlike the Lorentzian case, in Euclidean space we can consider closed 
worldlines, and actually these will be the ones relevant to pair production. A 
closed wold-line divides spacetime into two regions. Following [|J we denote 
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them as the "inside" and the "outside" of the worldline (see Fig. 1). Note 
that in flat space there is a natural way to assign these labels, but in de Sitter 
space (that is on the 2-sphere) these denominations are just conventional. 
We denote by "outside" the region in which, upon analytic continuation, 
the electric field does not change when the pair is created, maintaining its 
initial value E . By Gauss' law, the field in the inside region is given by 
Ei = E ± e. The double sign reflects the fact that the worldline has to be 
assigned a direction in which the charge flows, and in principle both directions 
are possible. Without loss of generality, we take the minus sign, 

E l = E - e, (8) 

since a change in the direction of the worldline amounts to a change in the 
sign of e. 

Using F^F^ = 2E 2 , we have 

2 

S E = M J s ds-(eE -j) J e. (9) 

Here V is the volume of the "inside" region. In we have dropped the 
constant term 

\ J Sx^glEl 

This term is just the Euclidean action for the background configuration, 
without the instanton. Of course this constant is physically meaningless 
because a constant can always be added to the action. It is customary 
to choose the constant so that the action for the background configuration 
vanishes [0. To summarize, the action @ is proportional to the length of the 
worldline minus a term proportional to the area enclosed by the worldline. 

2.2 Antisymmetric tensor field coupled to a membrane 

Let us now consider the natural generalization of (|3]), which describes an 
antisymmetric tensor field of rank iV 

A = A^___ p] dx^ A ... A dx p 

interacting with the N dimensional worldsheet of a charged membrane in a 
spacetime of dimension d — N+l [|J. To keep the discussion simple, we shall 
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consider the situation in which the field A is external, so that backreaction of 
the membrane on the field is ignored. This will facilitate also the comparison 
with the method of Bogolubov coefficients in 1+1 space-time dimensions, 
since with that method the electric field has to be treated as external anyway. 
The action is now given by 

S = -M Jy/=jd N £ + e J^A. (10) 

The first term is the Nambu action, proportional to the area of the worldsheet 
S, where 7 is the determinant of the worldsheet metric 7 a b, and £ a (a = 
0, N — 1) is a set of coordinates on S . In ( [TDD M is just a constant. For 
7V=1 this constant is the particle mass, for N=2 it is the tension of a string, 
and for N = 3 it is the surface tension of a membrane. The second term in 
( |T0D is the generalization of the electromagnetic coupling e / A^dx^ that we 
used in (§). Eq. fllCf ) can be Euclideanized using prescriptions analogous to 
the ones that led from (0) to (0), yielding 

S E = Mj^d N i + eJ^A. (11) 

As with the 1+1 dimensional case, the field strength associated with A, 
F = dA, has only one independent component, and can be written in the 
form (H), where now the antisymmetric symbol e has N + 1 indices. 

For closed worldsheets we can use Stokes' theorem and (|j) in ( p+| ) to find 

Se = M^ y^d N £ - eE J e, (12) 

where we have used a constant external electric field E Q . Without loss of 
generality we take 

E > 0. 

The integral in the last term is just the volume of the space-time region 
"inside" the closed worldsheet. Eq. ( |12"D has the same form as (§) if the term 



proportional to e 2 is neglected. Therefore one can neglect backreaction when 
|e| << Eq. 

Notice that (fl2]) is proportional to the area of the worldsheet minus a term 
proportional to the volume enclosed by the worldsheet. This has exactly the 
same form as the Euclidean action for the process of false vacuum decay 
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through bubble nucleation [fj] , in the limit in which the thickness of the wall 
separating the true from the false vacua is much smaller than the radius of the 
bubble. Both processes are similar in many respects [|J, the main difference 
being that membrane production by an antisymmetric tensor field can occur 
repeatedly at any given point in space, whereas vacuum decay occurs only 
once. 



2.3 The instantons 

The equation of motion following from (^) has been given for instance in 
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K = Y b K ab = (13) 

where K ab is the extrinsic curvature of the Euclidean worldsheet 

K ab = -e 6At V a n>\ (14) 

Here n M is the normal to the worldsheet, and = dbX^{E, c ) are the tangent 
vectors (our sign convention is that points towards the outside region.) 

In 1+1 dimensional flat spacetime, the only solution of (0) is a circular 
worldline of radius R = A4/eE , 

(4) 2 + 0* 1 ) 2 = Rl (is) 

Since we have chosen Eq > 0, in order for Rq to be positive we need 

e > 0. 

From (H), this means that the electric field inside the circle will be smaller 
than outside. 

The evolution of the pair after nucleation is given by the analytic contin- 
uation of (|1|) to Minkowski space 

- (x ) 2 + (x 1 ) 2 = Rl (16) 

This hyperbola has two branches. The branch on the right represents a 
particle of charge e moving forward in time. The one on the left represents a 
particle of charge e moving backward in time. This is interpreted in the usual 
way as an antiparticle of charge — e moving forward in time. The particle and 
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antiparticle pair nucleate at time x° = 0, separated by a distance 2R$ and 
with zero velocity. After that, due to the constant force exerted by the field, 
they start moving away from each other with constant proper acceleration 
Rq 1 . 

In higher dimensions there are also iV— spherical worldsheets which are 
extrema of the action. These represent the nucleation of spherical membranes 
[§]]. Let us consider directly the instantons in de Sitter space of radius H l . 
The flat space instantons can be obtained as the limiting case H — > 0. 

With a spherical ansatz for the worldsheet the action (|12"D takes the form 

S E = MS N (R ) - eE V N (6 ). (17) 

Here 

S«{«o) = ^TT< (18) 



is the surface of a worldsheet of radius Ro, 9o is the polar angle on the 
d-sphere of radius H -1 (see Fig. 1) such that 

R = H' 1 sin 9 , 

and Vat(#o) is the volume of the d-sphere that is enclosed by the worldsheet 
of radius Rq: 

V N (9 ) = -^r#- (JV+1) T sin^ 9d6. (19) 

r (e±ij Jo 

Extremizing fllTD with respect to 8q 

dS E _ 
d6 ~ 

we find 

tan# = iVtf^r, ( 20 ) 
eEo 

which means that the radius of the Euclidean worldsheet is 

r NM r 2 n 

(N 2 H 2 M 2 + e 2 E^y/ 2 ' y ' 
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Substituting ( BTj ) back into ( I7[) we find the Euclidean action for the instan- 
tons 



5, 



2nH 



-2 



{M 2 H 2 + e 2 El) l l % -eE 



AttH 



-3 



e£ 2#.M 
arctan 



S, 



2 e£ 



- 2eE 



(N=l) 
(N = 2) 



(N = 3) 



(22) 
(23) 



(24) 



_(9H 2 M 2 + e 2 E$y/ 2 

An interesting feature of Equations (22- |24]) is that one finds instantons of 
finite action both for e > and e < 0. 

For e > the electric field in the inside region decreases with respect to 
the initial value [see©]. From ( |2~0D this corresponds to 8 < |. For e < 
the electric field in the inside region actually increases with respect to the 
initial value. This corresponds to #o > f • m this case the "inside" region 
is actually larger than the "outside" one (see Fig. 1). Strictly speaking, 
if the electric field is treated as external, we should not say that the field 
increases or decreases in the inside region. However, both cases should still 
be distinguished. For instance, in 1+1 dimensions, the instanton with e > 
corresponds to the creation of a pair with the "screening" orientation. That 
is, after nucleation, the + charge is to the right of the inside region and the 
— charge is to the left (recall our convention E > 0). On the other hand, 
for e < the pair has the "anti-screening" orientation, with the + charge 
to the left and the — charge to the right. Similarly membranes can nucleate 
with two different orientations depending on the sign we take for e. 

It might appear that the particles in pairs with the anti-screening ori- 
entation would move towards each other after nucleation, and eventually 
anihilate each other. However, as we shall see, the distance between both 
particles actually grows with time due to the inflationary expansion. We 
should emphasize that the anti-screening instantons are just as physical as 
the screening ones. In order to find agreement with the results obtained 
using Bogolubov transformations, both instantons will have to be included. 
In the context of vacuum decay, the case e > corresponds to the ordinary 
transition from false to true vacuum, whereas the case e < corresponds to 
the decay of the true vacuum through nucleation of false vacuum bubbles [|] . 

In the limit when the electric field is switched off, Eq — > 0, the instantons 
become spheres of maximal radius Ro — > H" 1 , and the action (22-24) reduces 
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to 

S E = MS N {H- 1 ). (25) 

These are the instantons for the spontaneous nucleation of defects during 
inflation ||, which we shall consider in more detail later on. 

In general, for finite Eq, the action for cEq > is always smaller than 
that for eEo < 0. This means that it is more probable to nucleate a screening 
membrane than an anti-screening one, in agreement with naive expectations. 

In the flat space limit H — > 0, the anti-screening process has infinite 
action, and therefore it is not possible. Only the action for the screening 
instanton e > remains finite. For e > and H — > we have 

and 

M 

This expression reproduces the thin wall instanton action for vacuum decay 
in flat space |2[], where Ai is the tension of the wall and eE is the difference 
in energy density between the false and true vacua. 

2.4 Analytic continuation 

The evolution of the membranes after nucleation is given by the analytic 
continuation of the instantons back to Lorentzian signature. We have seen 
that, for the d — 1 + 1 case in flat space the instanton is a circle, and the 
analytic continuation is a hyperbola representing the world-line of the particle 
and antiparticle accelerating away from each other. Note that the hyperbola 
(|16D is centered at the origin x l = 0, but of course pairs can nucleate at other 



locations too. If we act on the instanton flT5|) with a space-time translation, 
the resulting trajectory 

is also an instanton. Thus (|15|) is just one solution out of a two parameter 
family. By analytically continuing and cle —> ia, we obtain a two 

parameter family of Lorentzian solutions 

- (a; - a) 2 + (x 1 - bf = Rl (28) 
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which represent pairs nucleating at any space-time point x^ = (a,b). 

Similar steps have to be taken to analytically continue the instantons 
describing the creation of membranes in de Sitter. As mentioned before 
these are N— spherical worldsheets of radius R embedded in d-sphere of 
radius H^ 1 . The instantons can be represented as the intersection of the 
d— sphere with a hyperplane at a distance 

UJq = H~ x COS^o 

from the origin, where 9q is given by (|20|), see Fig. 1. In the representation 
(H), the instanton is given by 

J=l 

X d = u , (29) 

where X A are the coordinates in a 'fictitious' embedding Euclidean space. 

Following H, to analytically continue this solution we choose the flat 
Friedman- Robertson- Walker (FRW) coordinates in de Sitter space, (t,x). In 
these, the metric takes the form 

= _ dt 2 + e 2Ht ds 2_ (30) 

These coordinates are related to X A through the relations 

X° = H- 1 sinh Ht + -Hx 2 e Ht , 

2 

X d = H' 1 cosh Ht - l -Hx 2 e H \ (31) 
X = xe m , 

where the vector X has components X J , (J = 1, d — 1). The coordinates 
(t, x) cover only half of the hyperboloid (0). 

Taking X d = ujq in ([H]), the worldsheet of the membrane after nucleation 
is given by 

x 2 = H- 2 (l + e~ 2Ht )-2H- 1 iu e- m . (32) 
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This solution represents a spherical membrane which is expanding in time, 
with physical radius given by 



R 2 = H -2 {e 2m + 2H - 1 u e m . (33) 

Note that sign(a;o)=sign(e), but R never vanishes for either sign of e. In both 
cases, the radius grows like the scale factor at late times. For N — 1 the 
spherical "membrane" reduces to a pair of points, whose worldline is given 
by x = ±(x 2 ) 1//2 , with x 2 given by (j32|). In this case R is one half of the 
physical distance between the particle and antiparticle in the pair. 

As with the flat space case discussed above, the solution (|32|) belongs to 
a family of solutions which can be obtained from a ci-parameter family of 
instantons. This family is obtained by applying 0(d + 1) rotations to the 
instanton (f2l^). The group 0(d + 1) has d(d + l)/2 generators. Of these, 
d(d — l)/2 leave the instanton invariant. They correspond to rotations in 
the space (X E ,X). The remaining d generators correspond to rotations in 
the (X d , X E ) plane or in any of the (X d , X 1 ) planes (i = 1, d— 1). These 
generators which do not leave the instanton invariant are the so called zero 
modes. Their effect is to rotate the hyperplane X d = u>o in (|29|) , effectively 
translating the center of the worldsheet to a new location on the rf-sphere. 

Upon analytic continuation the parameters corresponding to rotations 
in the (X E , X J ) plane (J 7^ 0) have to be complexified along with X E in 
order for the resulting solutions to be real. Recall that even in flat space, 
the parameter had to be complexified to obtain (1231). In the present 



case, rotations turn into boosts in the (X°,X J ) plane when the angle of 
rotation is complexified, 

ctE = iot. (34) 

In this way, the group of rotations 0(d + 1) becomes the group of de Sit- 
ter transformations 0(d, 1), which can also be thought of as the group of 
Lorentz transformations in the Minkowski space in which the hyperboloid 
(H) is embedded. 

Let us consider the case d — 1 + 1 in some detail. The general worldline 
after nucleation is obtained by taking X E — > iX° in (p9|) and then applying 
a boost in the (X°,X 2 ) plane followed by a rotation in the (X l ,X 2 ) plane 

X'° = X° cosh a + X 2 sinh a (35) 

X' 1 = -(X°sinha + X 2 cosha) sin (3 + X 1 cos/3 
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X' 2 = (X° sinh a + X 2 cosh a) cos (3 + X 1 sin (3. 

Here a and /3 are arbitrary parameters. 

Elliminating X 1 and X° from the previous equations one has 

X' 2 cosh a cos (3 — X' 1 cosh a sin /3 — X sinh a = X 2 . 

Taking X 2 = ujq, dropping primes and using the transformations fl3~T|) we 
find, after some algebra, that the general world-line after nucleation is given 

by ' 

R 2 = e 2^ (x _ Xq) 2 = H -2 {e 2H(t~t Q) + 1} _ 2#-i Woe JT ('-*o) > (36) 



where 



x = jj-i cosh a sin/? 

cosh a cos /3 + sinh a 



t = H ln(coshacos/5 + sinh a). (38) 

Eq. (|36|) represents a pair centered at the point xq. The parameter t shall 
be refered to as the time of nucleation. For Rq ~ H^ 1 this is somewhat 
conventional, since there is no precise instant of time at which the pair nu- 
cleates [H]. The problem is that the concept of simultaneity becomes blurry 
at distances comparable to the horizon. Strictly speaking, all we can say is 
that solutions with different values of to are time translations of one another. 
From a geometric point of view, however, it is clear that (xo,to) represents 
the center of symmetry of the worldsheet, and in this sense it is natural to 
think of it as the nucleation event. 

Note also that there is no absolute value sign in the logarithm in (|58|) . For 
u 7^ 0, the solutions with cosh a cos f3 + sin a > are qualitatively different 
from the ones with cosh a cos (3 + sinh a < 0. Actually, we shall see that the 
latter ones are unphysical. They correspond to pairs whose "inside" region is 
centered at spatial infinity, so upon nucleation the electric field would change 
over an infinite (and disconnected) region of space. To find agreement with 
the Bogolubov method these solutions have to be discarded. 

The same arguments can be repeated for N > 1. The general Lorentzian 
solution is a spherical membrane of physical radius (|33|) centered at any 
spacetime point B. 
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3 Perturbations and zero modes 



To compute the contribution of the instantons to the partition function we 
need to study small fluctuations of the instanton worldsheet. For this it is 
very useful to adopt the covariant formalism developed in Refs. [fill [TTL O, 



14j| , according to which the worldsheet perturbations are represented by a 
scalar field "living" on E, which has the meaning of normal displacement. 

Denoting by the coordinates in de Sitter space, the instanton config- 
uration will be denoted by x M (£ a ). Consider now a slight deformation of the 
worldsheet 

^(C)=^(C) + ^(C). (39) 

Since only deformations orthogonal to the worldsheet are physically mean- 
ingful, we can set 

fcc"(f a ) = .M- 1/2 0(rK. (40) 

Here n M is the normal to the worldsheet and 0(£ a ) has the meaning of a 
normal displacement. The factor Ai^ 1 ^ 2 is inserted so that (f) has the correct 
dimensions for a scalar field in N dimensions, [0] = (mass)^ 2 - 1-1 . 

Actually, the equation of motion for <fi can be derived from kinematical 
considerations. Since is a scalar on the worldsheet, it has to satisfy a 
covariant equation. The only tensors available in E are the metric 7 a & and 
the extrinsic curvature K a b, but because of the symmetries of our problem 
they are proportional to each other K a b oc 7 a f, (see e.g. [|[). The only 
covariant second order differential equation that we can write down with 
such ingredients is 

- A(f) + M 2 (f) = 0, (41) 

where A is the Laplacian on the spherical worldsheet. By symmetry, M has 
to be constant. 

To determine the value of M we can use "known" solutions of fl4"T|): the 
zero modes. These are field modes for <ft which do not correspond to true 
perturbations of E, but to the infinitessimal version of the rotations consid- 
ered in the previous section. These solutions can be found from geometric 
considerations. Let the instanton be given by [see (p9|)1 

(X° E ) 2 + ir(X J ) 2 = H- 2 (42) 
j=i 
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where, as usual, X A are cartesian coordinates in the flat space in which the 
Euclidean de Sitter space is embedded. 

The vector n M , orthogonal to S, can be thought of as a vector in this 
embedding flat space, with components n A . As such, it is tangent to the 
d— sphere of radius if -1 centered at the origin, and orthogonal to the N— sphere 
of radius Ro centered at X d = ujq, X 1 — (J ^ d). It is easy to see that the 
normal vector is given by 



n 



HR (Xloo, — Rl). 



Here X has components X 1 , (J ^ d). A small rotation of angle a in the 
(X J , X d ) plane induces the change 

5X J = aX d , 

SX d = -aX J . 

This transforms £ into a new worldsheet which is also a solution of the 
equations of motion. Therefore, taking X d = u in the equations above, the 
field 

AT 1/2 0(C) = n A SX A (C) = aiHRoy'X^C) (43) 

has to be a solution of (fill) , fo r any J = 0,...,d — 1. As is well known, 
the cartesian components X J (^ a ) of the points on the N— sphere are linear 
combinations of the spherical harmonics with L — 1. The spherical harmonics 
are eigenfunctions of the Laplacian with eigenvalue 

X L = -L{L + N - l)i?o 2 , (L — 0, oo) (44) 



so, taking L = 1, we have A0 = —NR 2 <p. Comparing with (fil]) we have the 
mass that we were looking for 

M 2 = —NRq 2 . (45) 

For eEo — > this reduces to M 2 = —NH 2 , a result which was found already 
in §. 

Of course, equations ([y]) and (|45|) can also be derived from a perturbative 
expansion of the action. Introducing (^) in the action and expanding to 
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second order in 0, the action for the perturbed worldsheet £ can be written 

as |g, m 0, 

S' E [E] = ^[E] + 4 2) [0], 



where 5b [E] is the action for the unperturbed instanton (|2^ - p4]) , and 



4 2) M = \j ^eV7(-0A0 + M 2 2 ) 



whith A the Laplacian on E and 



M 



Tl {N) - R$(gv - n»n v ) 



eEo 
M 



(46) 



(47) 



Here TZ^ N ' is the Ricci scalar on E, and R^J and g^ u are the Ricci tensor and 
the metric in the embedding de Sitter space [actually Eqs. 



and (|47]) are 

valid for pertubations to any worldsheet solution, embedded in an arbitrary 
curved spacetime of dimension d = N + 1]. In de Sitter space (i.e. on the 
c?-sphere) R^J = H 2 (d — l)g^ u , whereas on the iV-spherical worldsheet of 

radius R , = i ab R^ = N(N - l)Ro 2 . Substituting in (0) and using 
(|2lD for R , the effective mass M 2 simplifies to 



M 



-NR^ 2 , 



in agreement with QI5|). 

For later convenience we expand an arbitrary perturbation <p in terms of 
the (real) spherical harmonics on the iV-sphere, 4>lj 



LJ 

These satisfy A(f> LJ = Xl4>lj, with X L given by ([ 



and 



(48) 



(49) 



The index J labels the degeneracy for given L 

J = Q,...,g L -l, 

where 

9l 



(2L + N -l)(N + L-2)\ 



L\{N -I) 



(50) 
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For L = 1 we have, with the normalization (EDI 



N+1 



1/2 



X J (0, (51) 



where Sjy is given by ([Of). 

Comparing (|43D with (E8| ) and using (|5lD we find that for an infinitessimal 
rotation of angle cfctj in the (X J , X n ) plane 



'MS N {R 



1/2 



dC u = H-Haj ^7 + V 7 j • (52) 

This equation is often refered to as the normalization of the zero modes, and it 
will be important in order to interpret certain divergences in the semiclassical 
evaluation of the partition function. In the flat space limit, the square root 
in the right hand side of fl52"l) reduces to the familiar expression S^ 2 @. 



4 Related instantons 

The formalism used in the previous sections can be easily extended to study 
a few more instantons which are closely related to the ones we have seen so 
far, and which will be relevant for future discussion. 



4.1 Nucleation of topological defects during inflation 

We saw in Subsection |2.3| that for H ^ the action for the instantons 
remains finite when the external field is switched off. This corresponds to 
the spontaneous nucleation of membranes (or topological defects), due to the 
gravitational field alone. 

When the external field Eq is zero, there is no need to restrict ourselves 
to co-dimension 1. Thus monopole pairs, strings and domain walls can spon- 
taneously nucleate in d- dimensional de Sitter space, with d > N. The corre- 
sponding instantons are found by intersecting the (i-sphere with the necessary 
number of hyperplanes through the origin: 

X A X A = H- 2 (53) 

A=0 
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X l = 0, (i = N + l,...,n), 

which gives N-spheres of radius H~ l . (here we use the lower case latin index 
i for later notational convenience). 

When studying small perturbations around the instanton for co-dimension 
larger than 1, we will have more scalar fields 'living' on the worldsheet, one 
for each normal direction. The generalization of (|40|) is 

d 

5x"(C a ) = M~ 1/2 ]T 4>^{i a )n^. (54) 

i=N+l 

Taking the normal vectors to be perpendicular to the hyperplanes X 1 = 0, 
the effective action for <f>® is [||, [H], 



S (2) M4E fd^i^i-A + M 2 )^], (55) 



d 

L 

i=N+l ■ 

where M 2 is still given by (45) with Rq = H" 1 , i.e. 



M 2 = —NH 2 . 

The fields can be expanded in terms of spherical harmonics on the iV— sphere, 
0(0 — J2 Cj*j4>lji and the normalization of the zero modes can be worked 
out in exactly the same way as before. For an infinitessimal rotation in the 
(X J , X 1 ) plane (where J = 0, N, % — N + 1, d) of angle do$ we have 

dCf] = g-'dcS" ( ^fP j n '■ (56) 

Rotations of the X J among themselves or of the X 1 among themselves leave 
the instanton invariant, so these will not correspond to zero modes. Then, 
the total number of zero modes is (N + l)(d — N). 

In general, of the (iV+ l)(d — N) zero modes, d will correspond to space- 
time translations and the remaining N(d — N — 1) correspond to rotations 
in the spatial orientation of the defect 0. For instance, monopole pairs in 
3+1 dimensional de Sitter space can nucleate with all possible orientations 
of the relative position, and loops of string can nucleate with all possible 
orientations of the plane of the loop. 
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4.2 Particle at finite temperature in flat space 

For a massive particle at finite temperature the Euclidean action is 

S E = M f\l+x{t E )yi 2 dt E , (57) 

where t E is the Euclidean time and x(t E ) is the Euclidean trajectory, which 
has to be periodic x(t E ) = x(t E + /3). Here (3 is the inverse of the tempera- 
ture. The instantons for this system are straight lines wrapping around the 
compact temporal dimension 

x{t E ) = x = const. (58) 

The second variation of the action is found by substituting x(t E ) = x + 
M~ 1/2 4>(t E ) in (|57p, 

4 2) = ~ f{$?dt E . (59) 
2 Jo 

The fluctuations can be thought of as a set of massless scalar fields in 0+1 

dimensions. From (|5^), satisfies = 0. The zero mode solution = const. 
amounts to a spatial translation of the original solution 
Expanding (p in "spherical harmonics" on the 1-sphere 

= r 1/2 ^O+^E j^LCOS 

we see that the zero mode is in the L = sector. An infinitessimal translation 
dx corresponds to 

dC = dxo^/Mp. (61) 

The instanton (j58|) and the previous equation will be used to derive the 
partition function for a gas of massive particles at finite temperature. 

4.3 Pair production in more dimensions 

In curved space, E = const is a solution of the homogeneous Maxwell's 
equations only in 1+1 dimensions. However, in flat space one can take E = 
const in arbitrary number of dimensions. For completeness, we shall briefly 



2irt E 



D L sin 



27rti 



(60) 
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consider the well known 3+1 dimensional case in flat space. The action for 
a particle in an external field is still given by 



S E = Mjds + ej A^dx", (62) 

where we are integrating over a closed worldline. If the electric field E is 
in the x direction E = (E,0,0), we can take A y = A z = 0. It is clear 
from the previous sections that there will be an instanton solution which is 
a circle of radius Ro = M.{eE)^ x in the (is, a;) plane, where t E is Euclidean 
time. To find the second variation of the action, we take a perturbation to 
the instanton of the form Q51|). One of the normal vectors is chosen 
in the outward radial direction from the origin of the (t E ,x) plane. The 
corresponding perturbation is denoted by cj) r (s), where s = Ro&rctaii(t E /x) 
parametrizes the worldline. The second variation of S E with respect to <p r 
is independent of the existence of the two orthogonal directions y and z, so 
from subsection a, the quadratic action for r is that of a scalar field of mass 
M 2 = —Rq 2 in 0+1 dimensions. 

The other two normals are chosen in the y and z directions. The as- 
sociated perturbations <fi y (s) and <f) z (s) do not change the second term in 
(|B2"D, since A y = A z = 0. Therefore, from the previous subsection these 
perturbations will be massless fields. Putting it all together 

2 Jo 



Ro K) +r y + r z ■ (63) 



The zero modes for the massless fields are in the L = sector. They corre- 
spond to translations in the y and z directions and can be normalized from 
(§T|), with (3 = 2ttR , 

(y) ah^) 



dCty dC q 



dy dz 



2nR M, (64) 



The zero modes for r are in the L = 1 sector and they can be normalized 
from (p2) 

where we have replaced H~ x daj by the corresponding flat space limit dt E or 
dx. 
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5 Semiclassical partition function 



For systems at finite temperature, the lifetime of a metastable state is related 
to the imaginary part of the free energy [0, |2|] [see Eq. (|86| ) below]. The free 
energy is defined as 

F = -p-HnZ, (66) 
where f3~ l is the temperature and Z is the partition function 

Z = tr[e-?% (67) 

with H the Hamiltonian of the system. The key to the semiclassical evalua- 
tion of Z is to first express it as a path integral. For instance, for the case 
of a single non-relativistic particle in flat space moving in a potential V(x), 
one has (see e.g. ref. for a nice discussion) 



Z= [ Vx(t)e- SE , (68) 

Jx(0)=x(P) 

where the integral is over all paths which are periodic in Euclidean time with 
periodicity /3, and 

Se= I dt[\x 2 {t)+V{x{t))\ : (69) 
J o 2 

is the Euclidean action. 

As is well known [T6L de Sitter space behaves in some respects like a 



system at finite temperature (3~ l = H/2tt. One difference with flat space 
is, however, that on the sphere all directions are compact. Therefore, once 
in Euclidean there is no actual distinction between temporal and spatial 



directions. In our case, the Euclidean action is given by (12) and the natural 
generalization of fl6"8"|) is 



PS(C)e" 5Ep] , (70) 

where now the integral is taken over all closed worldsheets (or worldlines for 
N = 1). 

In the semiclassical limit, Z will be dominated by the stationary points of 
Se, and so it will be a sum of contributions from multi-instanton configura- 
tions Z = Zq + Zi + ... . Here Z^ is the contribution of a configuration with k 
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widely separated instantons. This configuration has action IcSe, where Se is 
the action for one instanton. Hence Zk will have the exponential dependence 

Z k oc e~ k§E . (71) 

To find the preexponential factor, one has to integrate over small fluctuations 
around the stationary points. Eq . (|70|) is rather formal, because we have not 
specified the measure of integration (this can actually be quite complicated, 
due to reparametrization invariance.) However, if all we are interested in 
are small fluctuations around the instanton solutions, the integration over 
neighboring worldsheets amounts to an ordinary path integral over the per- 
turbation fields 0(£ a ) that we introduced in section [3]. Therefore, we can 
write 

oo e -kS E ( ( _ _ s(2) 



Z = f —{J V(Pe ) ' (?2) 

where S E is the second variation of the action, given in (|46|). The sum 
is over multi-instanton configurations. The path integrals over give the 
contribution of fluctuations around each one of the spherical worldsheets. 
The k\ in the denominator can be understood as follows ||. When integrating 
over we are integrating also over the zero modes. That means that we are 
integrating over all possible locations of the instantons in Euclidean space. 
Since the instantons are identical, we must divide by k\ to avoid overcounting. 
Eq. (|72|) can be rewritten as 



Z = e z \ (73) 

where 

Z x = e- §E J V<P exp (y 0(A + M 2 )<P^d N ^j . (74) 

This is just the path integral for a free scalar field on a curved background 
(the sphere), and can be calculated using well known manipulations. 

Following flT7| , the field is expanded in spherical harmonics (fiij [see (|4~8D]. 
The integral over can then be expressed in terms of the coefficients Clj 

^=n^- (75) 
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Note that Z\ is dimensionless, whereas Clj has dimensions of (mass)^ 1 . To 
render T>ip dimensionless one has to introduce the parameter fx with dimen- 
sions of mass. Using (HSI) and one has 




Z x = e~ s * J] / ^777^ exp --[£(M 2 - A L )CL] , (76) 



where A^ is given by flU]) . After Gaussian integration one obtains 

Z t = e"*- I][(^o)AZ 1/2 ] 9L = e- s -(det[(Mflo)- a O])- 1/a , (77) 

L 

where gr^ is given by (0). Here we have introduced the dimensionless oper- 
ator 6 = Rl(-A + M 2 ), with eigenvalues A L = Rl(M 2 - A L ), [see (|§]. 

Since the eigenvalues are known, the determinant can be calculated with 
the ^-function regularization method. In terms of the generalized ^-function, 

az)=Y.9LA- L z , (78) 

L 

the determinant can be expressed as 

detK/i^or'O] = QiRoyXWe-t'W. (79) 

The values of ((0) and ('(0) for the operator O on the iV-sphere are cal- 
culated in the Appendix. For N = 1 we obtain £(0) = and C'(0) = 
-21n[2sinh(7i\R M)], so that 

p— Se 

Z x = — r. (N = 1) (80) 

Note that, since ((0) = 0, the dependence on the arbitrary renormalization 
scale fx disappears. 

For worldsheets of dimension AT > 1 it turns out that ((0) = only for 
odd N \TE\. In general, for even N the determinant will depend on the renor- 



malization scale fx. As mentioned before, the calculation of Z\ is the same 
as the calculation of the effective action for a free scalar field in a curved 
spacetime of dimension N. The appearance of a renormalization scale in 



that context is well known [1^, O] (i n particular, this scale is responsible 
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for the so-called trace anomaly, the quantum mechanical breakdown of con- 
formal invariance). The difference between even and odd dimensions can 
be understood from the fact that in dimensional regularization the infinities 
come from the poles of the gamma function T[j — (N/2)], where j is an in- 
teger. As a result, for odd iV the effective action is finite after dimensional 
regularization, and the usual log \i terms do not appear. In general, for even 
iV we have to live with the arbitrary scale /i, unless ((0) happens to vanish 
accidentally for our particular worldsheet geometry and mass of the scalar 
field <p. It is well known that ((0) can be computed in terms of geometrical 
invariants. For N = 2 we have 



(81) 



where m 2 + £1Z^ is the quantity that we have denoted as M 2 [see Eq. ([47])]. 
Notice that both in flat and de Sitter backgrounds, the term involving the 
external Ricci tensor in ([Tj]) is constant, and can be included in m? . We shall 
come back to the discussion of /x and its role in the context of renormalization 
in the next Section. 

Our derivation of the semiclassical partition function has been rather 
formal. As a check, let us apply these ideas to a simple example where the 
result is known: the case of free particles at finite temperature in flat space. 
This will also illustrate the general procedure for dealing with the zero modes. 

For particles at finite temperature, the action for the instantons (discussed 
in subsection |4~2] ) is equal to the mass M. of the particle times the length of 
a worldline that wraps around the compact temporal dimension Se = M.f3. 
The perturbation field <fi is massless, M = 0, and so the expression ( |30D di- 
verges because of the vanishing denominator. This is to be expected because 
the operator O has a zero eigenvalue corresponding to the translational zero 
mode discussed in section |3] (for simplicity we start by considering only one 
spatial dimension transverse to the worldline, hence only one field <fi and one 
zero mode). 

Noting that (2n)- 1 ' 2 J dC exp(-M 2 C%/2) = M" 1 , the divergence at 
M — > can be avoided by leaving the integral over dCo undone and by 
excluding the factor M _1 from the determinantal factor. Then we can write 

dZ, = ^lH)- 2 6])-^ e -^^, (82) 
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(2vr.Ro) 2 (83) 



This can be cast in a more familiar form by setting T = f3~ l and increasing 
the number of transverse dimensions to three. Each transverse dimension 
brings an additional power to the preexponential factor. Interpreting d 3 xo 
as the volume element we have 

This is the correct expression for the "one particle" partition function of an 
ideal gas with the Maxwell-Boltzmann distribution [note that the instanton 
method is only valid when the exponent in ( [34]) is large, in which case there 
is no difference between bosonic and fermionic distributions]. The grand 
canonical partition function is obtained, according to (|73|), by exponentiating 
this expression. 

Actually, Eq. ( [73] ) gives the partition function for the case of vanishing 
chemical potential, something that we have tacitly assumed in our derivation. 
The effect of a chemical potential jl is to replace Z\ in (|73|) by e^Z\. The 
number of particles M in the volume V is then given by 

For vanishing jl, we have 

M=Z l . (85) 
Therefore Z\ is also the equilibrium number of particles. 



where 

det'^rd] es lim de *[W" 2 6] = 

is the determinant without the zero eigenvalue. 
Using fl5Ip and 2nR = [3 we have 
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6 Nucleation rates in flat space 



In flat space and at sufficiently low temperatures, the decay rate of a metastable 
state is given by [|7| 

T = 2\Im F\ = 2f3~ 1 \Im Z t \. (86) 

(For temperatures [3~ l > Rq 1 , this formula has to be modified 0. We defer 
the consideration of the high temperature regime to a later paper.) As shown 
in the previous section, the calculation of Z\ reduces to the calculation of 
a functional determinant. The fact that F has an imaginary part is due to 
the fact that the action in ([K]) has a negative mode, corresponding to 
L = 0, so the determinant will be negative. Upon taking the square root a 
factor of % will emerge. 

The membrane creation rate can be expressed as 

dT = ~ x 2p- 1 \det'[(fxR y 2 6]\- 1/2 e- §E \J\dVdt E , (87) 
where the Jacobian is given by 

|J| " dVdT E • (88) 

This equation follows from (|8~6"D and (|77|) . As before [see (|83|) 1 the prime in 
the determinant means that the (N + 1) zero modes (which are now in the 
L = 1 sector) are omitted, because the integration over dC\j is left undone 

The Jacobian, which is needed to change variables from dC\j to the Euclidean 
space-time volume element dV dt E , can be read off from (|5^| ) [for obvious 
reasons, in the limit H — > we replace H~ x daj by dt E (for J = 0) or by dx 
(for J = 1)], 



JV+l 

( MS N (R ) \~~ 



The overall factor of 1/2 in the r.h.s. of (|87|) is explained in Refs. |20], g]. 
It arises because the free energy of an unstable state can only be defined by 
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analytic continuation from a Hamiltonian in which the same state is stable. 
As a result, the contour of integration over the negative mode dCo has to be 
deformed into the complex plane in such a way that only half of the saddle 
point contributes to Im F . We should note, however, that in the derivation 
of (jS6|) given in [Tj, these considerations are not really relevant; and 2\Im F\ 
is essentially a convention to denote the r.h.s. of (|87|). 

Integration over Euclidean time cancels the factor of and we are 
left with the rate per unit volume 



dU r (MS N (R ) 



JV+1 



dtdV V \2tt(N+1) 



\det'[(fiR )- 2 O}\- 1/2 e- bE , (90) 



where Af is the number of membranes created. Let us now evaluate this rate 
for space-times of different dimensionalities. 



6.1 Pair creation in 1+1 dimensions 

From (P5| ) and ( |A7| ) the primed determinant is 

4sinh 2 (7ri? M) 



M 2__ R 2 (M 2 + i?o 2 ) 2 



71 Hq, 



Se 



hence, from (|90|) the creation rate per unit length is 

r m 

L ~ 2ttR ( 

Taking R Q = A4/eE , we have 

T eE ( 



2tt 



exp 



(91) 



This can be compared with the results of Stone f2lfl , who computed the 
vacuum decay rate in the in the sine-Gordon theory with non-degenerate 
vacua. His one loop result was 



r 

L 



eEc 
2vr 



\ln(l 



_ttM± 

e eE o 
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where M. is the mass of the kink and eEo is the vacuum energy density gap 
between neighboring vacua. As expected, the instanton calculation gives a 
good approximation when the Euclidean action is large Se » 1- 

Eq. (|91]) also gives the rate at which a metastable cosmic string will break 
up by nucleating pairs of monopoles || . In that case eE should be replaced 
by the string tension and M. by the mass of the monopoles. 



6.2 Pair creation in 3+1 dimensions 



As explained in subsection [4.3| , in 3+1 dimensions in addition to the radial 
perturbations of mass M 2 = — -Rjj, we have the perturbations <p y and <p z which 
are transverse to the plane of the instanton. These behave like massless fields 
M 2 = 0. Each field contributes its own determinantal factor, which for cf) y 
and (j) z is given by (|83|). Also, from (p4]), each contributes a factor (RqM) 1 / 2 
to the Jacobian. With this, Eq. ([91] ) is modified into 

T (eE ) 2 -?mL , \ 

v = - • < 92) 

This coincides with the rate of production of charged bosons in scalar elec- 
trodynamics 

in the limit Se » 1- (For a more thorough account of monopole production 
by a magnetic field, and pair production in the strong coupling regime see 



Refs. [22, E3 



6.3 String creation in 2+1 dimensions 

From @g) and (TAT0|) , 



^ = (M^y\^ R -s e -s E _ (94) 

Note that, since ^(0) ^ 0, the determinant depends explicitly on the renor- 
malization scale /i [the exponent 7/3 can also be derived from (|81|).1 Because 
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of the arbitrariness in fi we cannot give an absolute estimate of the nucleation 
rate the way we do for N=l or N=3. However, since fi is a constant, we can 
still predict how the rate changes when we change the external field E Q (that 
is, when we change i?o-) It is seen that the dependence of the prefactor on 
Ro is more complicated than what one would have guessed form dimensional 
analysis. 

So far, by using the ^-function method, we have avoided the question 
of ultraviolet divergences, since they are automatically removed by analytic 
continuation [[19] (see also Ref.|24| for a recent discussion, and references 



therein). However, we should recall that the functional determinants contain 
such divergences, and that these can be elliminated by suitable counterterms 
in the action. For N = 2, all divergences can be removed by counterterms of 



the form |19| 



c J d 2 ^+dj d 2 ^R. 



The first term is a contribution to the membrane tension. The second is a 
topological invariant which does not contribute to the equations of motion. 
Note, form (^) and (|8T|), that a rescaling of the arbitrary parameter [i can 
be reabsorbed in a redefinition of c and d. Then, we can elliminate the 
renormalized c by rescaling fi, and the renormalized d by shifting the string 
tension A4. 

In Ref. a different approach was followed in which the product over all 
eigenvalues was cut off at some physical scale /i. This method also produced 
the factor (^Rq) 7 ^ 3 , which was refered to as the "universal term" |25|, pEfl . 
The introduction of a physical cut-off is not unreasonable when dealing with 
an effective theory. In deriving the action for the perturbations we have only 
expanded to second order in (f>. This is justified at low momenta L, but in 
the ultraviolet limit, higher order terms become important. To realize that 
this is so, let us take the case of a flat membrane and Eq = (this is actually 
quite representative of the general case). The exact Nambu action is then 



S E = M j d^^fl + M-^da^). 



When d a (p is of order .M 1//2 higher order terms become important and the 
quadratic approximation fails. This suggests taking /i ~ M 1 ! 2 as a cut-off, 



which is essentially what was done in [25 
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6.4 Membrane creation in 3+1 dimensions 

Using ([50]) and ( A14 ) we have 



I = (gW)'*SU<-V*. (95) 

This can be compared with the results of Affleck Pfifl . He studied the decay 
rate of false vacuum in the theory of a scalar field with a symmetry breaking 
potential and nondegenerate vacua, in the limit in which the thickness of 
the wall separating the true from the false vacuum is much smaller than the 
radius of the bubble at nucleation. 

The first factor in the r.h.s. of (|95|) is the Jacobian |J| that comes from 
the normalization of the zero modes. The rest of the pre-exponential factor 
coincides with what Affleck calls the "universal terms". These are due to 
fluctuations of the worldsheet and therefore are independent of the details 
of the field theoretic model. That is, they are independent of the internal 
structure of the "kink" or domain wall. The evaluaton of quantum correc- 
tions to the effective action due to finite thickness of the wall is in itself an 
interesting subject, and the corrections can be important in realistic theories 
[p6| , p2| |23"1 , |27fl . However, since these corrections are not essential to the 
physics of nucleation here we shall concentrate on infinitely thin membranes, 
without internal structure. We will only mention that the fluctuations around 
thick walls can be treated very naturally within the formalism described in 
this paper. In addition to the tachyonic field of perturbations that we have 
considered, there will be an infinite tower of massive fields <p q 'living' in the 



N-sphere p8| , |TT|| , each one contributing its own determinantal prefactor. We 
plan to return to this question in a later paper where the thermal properties 
of the fields <p q (with temperature 27ri?o) are emphasized. |29l 



7 Nucleation rates in curved space 

Although the calculation of determinants for the instantons in de Sitter space 
offers no special problems (we only need to substitute the appropriate values 
of Ro in the expressions found in the Appendix); the definition of a nucleation 
rate in de Sitter is a more subtle issue which has often been eluded in the 
literature. One possibility would be to simply use Eq. (p0|), interpreting dV 
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as the physical volume element at the time of nucleation dVo 



dT = \X\dV , (96) 

where |A| is defined as the r.h.s. of (|9"Q). 

Although we believe that this expression is correct (when properly inter- 
preted) it clearly needs further justification. First of all, the physical volume 
element, dVo, is proportional to a power of the scale factor e Ht ° at the time 
of nucleation. If Eq. (|90|) was derived from a purely Euclidean calculation, 
how does the exponential of a Lorentzian time find its way into the r.h.s. of 
(IHD? Also, as pointed out in f|, the time of nucleation in de Sitter space is a 
somewhat ambiguous concept when the size of the instantons is comparable 
to the horizon, and in principle it is not clear what time one should use in 
dV . 

These difficulties prompt us to search for an alternative way of calcu- 
lating nucleation rates in de Sitter space, which does not take flSUP as the 
starting point. Recalling that de Sitter space behaves in some respects like 
a thermodynamical system, one can try to estimate directly the equilibrium 
distribution of membranes. For this one can use Eq. fl85|) , with A/~ the number 
of membranes and Z\ given by (|74|) . It is clear that 

d-i 

dU= (det'[(fiR )- 2 O])- 1/2 e- SE H(2n)- 1 / 2 dC 1J . 

j=o 

Using (|52"D, we have 

d-l 

dN = \H- d [] da j 
j=o 

where 

\ N+1 

x ^{^wfi) ' ^Wd])- 1 '^. (97) 

Eq. ([52]) is valid only for infinitesimal rotations, and in that case Yldotj 
can be identified with the differential solid angle on the d-sphere, dfl, within 
which the center of the instanton worldsheet is to be found, 

dN = XH- d dn. (98) 

One can interpret this equation as the "equilibrium distribution of instan- 
tons" in Euclidean space. 
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Of course, an equilibrium distribution of instantons is not a measurable 
object. However, it can be given meaning by conjecturing that the equi- 
librium distribution of membranes in the Lorentzian section is given by the 
analytic continuation of the previous object to real time. This prescription 
is just heuristic, and we do not know how to justify it further, except by 
saying that it reduces to ( p0|) in flat space and that it is quite natural from 
the mathematical point of view. 

To see how the analytic continuation is done, let us consider, for sim- 
plicity, the 1+1 dimensional case (higher dimensional cases are completely 
analogous). Then (|98|) reads dM = XH~ 2 cos a Edasdp, where otE and j3 
are polar and azimuthal angles on the 2-sphere. Upon analytic continuation 
oie = ia [see fl5%D], 

dM = \\\H~ 2 cosh adad/3. 

Note that the factor of i from daE cancels the imaginary factor from the 
square root of the determinant in A, so that dM is actually real. Using Eqs. 
( [37D and (j38|), one can change variables from (a,/3) to (xo,to). The Jacobian 
is 

d(to,x ) 2 _ m 

H coshae °, 



d(a,(3) 
and therefore 

dM = \X\e Ht °dx dt . (99) 
Generalizing to spacetimes of arbitrary dimension we have 

dM =\X\e NHto dx dt 0} (100) 



which is very similar to (p6|) , but now all ambiguities in dV have been re- 
solved. 

An equation of the same form as the previous one was given in ||, based 
on kinematical considerations. However, the parameter A was left unspeci- 
fied. Like in ||, Eq. (|100|) is a distribution in the space of parameters xq 



and to, and it is therefore independent of the time of observation. 

Following , we can use ( |36| ) to express t Q in terms of the physical radius 
R, and thus find the size distribution of membranes (or bubbles), 

dM |A| R(R 2 - RlY l/2 JD 

jdR. 



dV phys tf^o + OR 2 -^) 1 ^] 
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where dV p h ys = exp(NHt)d N xo- Notice that the distribution diverges at the 
lower end R — > H^ 1 when u>o < 0. This divergence was interpreted in [||]. 
For large radii, one finds the scale invariant distribution 

dM IAI dR 



dV rby , H d R d 

which depends on the external field Eq and membrane tension only through 
the coefficient A. The fact that the size distributions are time independent 
is easy to understand. As the bubbles are created, they are stretched and 
diluted by the inflationary expansion, giving rise to a stationary distribution 
of sizes [^|. 

The "nucleation rates" |A| for d — 2,3 and 4 can be read off from the 
r.h.s. of (|9lD , fl94j) and fl95f ) respectively, where Rq is given by (|2lD and 
Se by ( |2"2| - p4"|) . This covers the case of co-dimension one, d — N + 1. For 
completeness, we shall also consider the case of strings and monopole pairs 
spontaneously nucleating in 3+1 dimensions. 

For the case of strings, the co-dimension is 2, and so there will be two 
independent perturbation fields and two determinantal prefactors. From 
fl94|) it is clear that A will contain one factor of Ai 3 and a factor of /i 14//3 , 
where fi is the renormalization scale. Taking Rq = H~ l , the rest will be a 
numerical factor (which can be absorbed in /z), times the appropriate power 
of H necessary to give A the dimensions of (mass) 4 : 

\\\=^M 3 H- 11 / 3 e-^. 

Because of the arbitrary renormalization scale, we cannot obtain an absolute 
estimate for |A|, but only its dependence on the expansion rate H. Like in 
subsection ( |6.3| ), to obtain a crude absolute estimate one can set /i ~ A4 1 / 2 , 
which gives 

>M\ 23/3 



|A| ~ (-1 H A exp(~ATrMH- 2 ). (102) 

For A4 » H this can be considerably larger than the naive dimensional es- 
timate |A| ~ .£f 4 exp(— Se) mentioned in ||, or even the more sophisticated 
|A| ~ M. 3 H exp(— Se)- Unfortunately, the existence of an arbitrary renor- 
malization scale leaves us quite uncertain as to the overall normalization of 
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Let us now consider the case of pairs spontaneously nucleating in 3 + 
1 dimensions. The co-dimension is 3 and so there will be 3 perturbation 
fields <f>( l \ each one of them with mass M 2 = —H~ 2 , contributing to the 
determinantal prefactor. Each field has 2 zero modes [in the L = 1 sector, 
see (1561)1) which makes a total of 6. Four of them correspond to space-time 
translations and two of them to changes in the orientation of the monopole 
pair in three dimensional space. Thus we have 



V 2tt J \H 2 

The factor of (J\AH/2it) 3 comes from taking the third power of the prefactor 
in (|9lD, with Rq = H -1 . The factor 4irH~ 2 is a correction due to the fact 
that two of the zero modes represent changes in the angular orientation of 
the pair, rather than translations. For each angular variable, the Jacobian 
fl88[) has an extra power of H in the denominator [see (|5^)]. Integration over 
all possible orientations gives the factor of Atc. To summarize, 

|A| = ^expf-^V (1=3+1) (103) 



2vr 2 "V H 

where M. is the mass of the monopole. 

For the case of pair creation a 'size' distribution such as ( |101| ) is not very 
useful. Instead, it is more convenient to find the momentum distribution of 
particles. Let us find the conserved momentum as a function of the time of 
nucleation t . In 1+1 de Sitter space, the vector potential 

= -H^E^b^ (104) 

represents a constant electric field (note that F fll/ F tiu = 2£q). With this, the 
action for the point particle coupled to reads 

S = -Mj(i- x 2 e 2Ht f' 2 dT - ^ J e m xdr, (105) 

where r is the proper time and a dot denotes derivative with respect to r. 
Since the Lagrangian does not depend on x, the momentum 

k = ^ = Mxe 2Ht - eE e Ht , 
ox 
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is conserved. Setting x = in ( |36| ) we have 

j e 



dx Huj n e~ H{t ^ - e~ 2Ht 



dt Hx 

from i 2 — x 2 e 2Ht = 1 

fore, 



and from i 2 -x 2 e 2Ht = 1, it follows that i = e m °H\x\(l - H 2 cu^y 1/2 . There- 



k = -^^e Hto sign(x), (106) 

where we have taken into account that the particle to the left of the inside 
region has charge of opposite sign. 

Using this equation in ( |99|) , with |A| given by the right hand side of (|9l|). 
we have 

dJv s^dk , , . 

= e s E d=l + l 107 

dx 2n 

In 3+1 dimensions (without electric field), using (|100D , ( |103|) and (|106|) , we 
lUV dM ( 2nM\ k 2 dk , , . . , 

Note that the momentum distributions are flat, as expected on the grounds 
of scale invariance. 

At any given time t we only need to integrate up to a cut-off momentum 

since particles with higher value of the coordinate momentum have not been 
created yet [see ( |106| )1. Then 

n~— ^-e~ §E , (d = l + l) 



and 



M 3 



where n is the number density of particles per unit physical volume. As 
noted in ||, for the case of vanishing electric field the distribution of particles 
contains a Boltzmann factor exp(— Ai/T) where T = H/2tt is the Gibbons- 
Hawking temperature [IS]. 
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8 Charged scalar field in 1 + 1 dimensional de 
Sitter space 

The Klein-Gordon equation for a charged scalar field <p coupled to an external 
electromagnetic field is 

- ^(V„ - ieAJ{V„ - ieA„)cp + M\ = 0. (109) 

As we shall see below, for a constant electric field in 1+1 de Sitter this 
equation can be solved in terms of special functions. Then the problem of 
particle creation is amenable for calculation using Bogolubov transformations 
(see e.g. [fL9lQ . 



To apply this method it is necessary to specify an 'in' state and an 'out' 
vacuum. The 'in' state, \in >, is the physical quantum state of our system, 
fixed by initial conditions. The 'out' vacuum, |0 > ou t, is a different quantum 
state in the Hilbert space, whose choice amounts to a definition of particles 
at late times. Whether or not one can unambiguously specify |0 > out depends 
on whether or not it is physically reasonable to define particles at late times. 
One way to guarantee a reasonable definition is to switch off the gravitational 
and the electric fields at late times (although this may not be necessary). 

To illustrate the procedure, consider the case of vanishing electric field 
first ||30|| . In terms of the conformal time rj = —H~ 1 e~ Ht , and with ip = 
ip k {j])e ikx , Eq. reads 

/ M 2 \ 

<& + [ k * + Jp-l) <P* = o. ( n °) 

where a prime denotes derivative with respect to rj. This equation is sym- 
metric in k, and following the usual convention, we take k > (For E = 0, 
the results for k < are the same). Eq. ( |110|) has the general solution 

MV) = (l) 1/2 [A k Hl 2 \k V ) + B k HW(k V )} , (111) 



where 

v= [\- m 



M 2 V /2 
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and if W are the Hankel functions. For the 'in' state, we shall take the Bunch- 



Davies vacuum [30]. This is characterized by positive frequency modes of the 
form (|m|) with A k = 1, B k = 



1/2 



H<?\kr,). 



(112) 



The choice of this vacuum as the physical 'in' state can be motivated from 
many different points of view, and it is a clear favourite in studies of inflation. 
In the open coordinate system that we are using, this is the only truly de 
Sitter invariant vacuum The two-point function in this state coincides 
with the Euclidean two-point function |RI|, which has the important property 



of having the Hadamard form (roughly speaking, this means that it has 
similar ultraviolet behaviour as the two-point function in flat space). Also, it 
is believed that if the Universe nucleated from 'nothing' into a de Sitter phase, 
then this is the quantum state that the fields would be in after nucleation 
§0. 

To define the 'out' vacuum it is convenient to write down the equation 
for the scalar field in terms of the cosmological time t, 



<p k + Hip k + [M 2 + 



a 2 (t) 



¥k = 0, 



(113) 



where a dot denotes d/dt and a(t) = exp(Ht). Introducing ip k 



-1/2 



have 



M 2 - — + 



4 



a' 



-<Pk 



0. 



4>k we 



(114) 



Let us take A4 » if and define t}. as the time when the physical wavelength 
a(t)/c _1 is equal to the particle's Compton wavelength 



ke 



-Hti 



M. 



(115) 



For t»% the k 2 a 2 term is negligible compared to M 2 . Then we will have 
approximate solutions of the form 



- 1/2 \a 



-iwt 



D k e 



\-iwt 



where 



w 



1 \ V2 

-H 2 

4 



(116) 
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Since A4 » H, the exponentials oscillate very fast compared to the rate at 
which a 1 / 2 changes, and so we will have an approximate definition of positive 
and negative frequency 'out' modes. For t » t k , 



(117) 



This definition is not very natural for Ai 2 << H 2 , so we shall not be con- 
sidering this limit. 

A purist would object that we cannot define particles unless the expansion 
of the Universe is switched off. Consider then a FRW model in which the 
expansion rate 7i = a/ a is time dependent. Then i/) k satisfies the equation 



,.2 K 2 k 2 1 • 
M ~- + ^--2 n 



0. 



118) 



After a sufficiently long period of inflation, a oc exp(Ht), the expansion is 
adiabatically switched off starting at time t s , in such a way that H << M 2 , 
until a(t) reaches a constant value. The space-time is Minkowski in the 
asymptotic future. It is clear from (|118 ) that if t s > t k , the mixing between 
positive and negative frequency modes will be negligible during the period in 
which the expansion is being switched off. This means that for modes such 
that t s > tk the number of particles that is calculated using the definition 
(|117|) for positive and negative frequency modes is the same that the number 
of particles that would be found in the 'out' Minkowski region. 

The 'in' positive frequency mode can be expressed as a linear combination 
of the 'out' positive and negative frequency modes 



(119) 



The coefficients a k and /3k are the so-called Bogoliubov coefficients. Using 
the asymptotic expression for the Hankel functions at late cosmological times 
t — > oo (i.e. ?7 — > 0) 



8J 



1/2 i 

Z/7T 




and the relation 



T(l - v)e- 

[i + oWl, 

K| 2 - |/3 fc | 2 = 1, 



5 L> 




T(l + u)e 



+1 



x 



(120) 
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one readily finds 

\p k \ 2 = [ex P (27iwH- l )-l]-\ (121) 

where w is given by Qll6| ). This equation was found e.g. in Refs. |35|, 
(see also references therein), using a different coordinate system. 

The number density of particles per unit coordinate volume is then 

d -^ = \P k \ 2 f- ^e-^^f , (d=l + l), (122) 
ax 2tt 2tt 

where we have used M. » H. The same calculation can be done in d 
spacetime dimensions. The only difference in the final result is that 



w = (m 2 - { -^h 2 



2 \ 1/2 



in Eq. ( |121|) . Using M 2 >> H 2 we have 

dU ,„ 2 d 3 k ( 2nM\ k 2 dk 



a* IAI (2^ Kexp l— J IP- (d=3+l) (123) 

Comparing ( |122| ) and ( |123| ) with ( |107| ) and ( IDS ) we find complete agreement, 
not only in the exponential behaviour, but also in the the preexponential 
factor. 

Let us now consider the case of non- vanishing electric field in 1+1 dimen- 
sions. In conformal time, the Klein Gordon equation with given by ( |104j ) 
reads (note that V ^ = 0) 

H 2 V 2 ^ + (k 2 HS 2 - 2eE k V + M 2 + ^f^j <p k = 0. (124) 

This is the Whittaker equation, having the general solution 

= A k W Ka (2ik V ) + B k W^ a (-2ik V ), (125) 
where W\ a are the Whittaker functions, with A = +ieE Q H~ 2 and 

U-w—h^) ■ (126) 
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At early times (77 — > —00), Eq. (|124 ) is very similar to ( |110| ), and both 
reduce to 

ip'l + k 2 ip k = 0. 

Using the asymptotic expression for Hjf>(kr}) at large rj, the 'in' positive fre- 
quency mode in the Bunch-Davies vacuum has the form (p^\ ~ (Ank)^ 1 ^ 2 exp( 
and so it is positive frequency with respect to the conformal time. With the 
electric field switched on, we would like to choose a positive frequency mode 
which has similar behaviour at 77 — > —00. Noting that 

and since W^ a {z) = W-\ z), it is clear that we have to take Bk = in 

(HID, 

<4tl = (^kr^W^kv), (127) 

where the normalization is due to the usual Wronskian condition. 

It is quite straightforward to show that the state defined by the modes 



( 127 ) is a Hadamard vacuum. For this, one simply writes the two point 
function as a sum over the modes ( |127| ). Similarly, one can write the Bunch- 
Davies two-point function as a sum over the modes (|112|) . Because of the 



similarity in the ultraviolet behavior of both sets of modes, it is easy to show 
that the difference between both two-point functions is finite in the limit of 
coincident points, which simply means that both two-point functions have 
the same singularity structure. 

To define particles at late times (r/ — > 0) we observe that, in this limit, 
Eqs. ( |1 10| ) and ( |124| ) are again very similar, so the procedure that worked 



there will work here too. In terms of cosmological time t and the variable 

ipk = a 1//2 <y9fc we have 

f kj2 e 2 E 2 H 2 k 2 2eE k\ , n 

Let us denote by the time at which the physical wavelength of the mode 
is equal to the effective Compton wavelength 

ke- Htk = Im 2 + -^\ . (129) 
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Provided that 



e 2 E? 



H 



2\ V2 



» H, 



(130) 



it is clear that for t » tk we can define particles in the mode k, using as 
positive and negative frequency modes the solutions 



(131) 



As before, these definitions are not meaningful for w << H, and we shall 
not consider this limit. 

Let us introduce the new Whittaker function M\ a , 



M Ka {z) = T(2a + l)e* 



-i7r(o-+|) _ 



WxJz) 



T(a + A + 



+ 



r(<7-A + D 



(132) 



3 1 

-7r < arc z < -n 
2 & 2 



For small z = 2ikrj we have [0 

M x>a = z^ +ff [l + 0(z)], 

and so this behaves like k in fll31| ) provided that we take a = +i\cr\ (recall 
that a is pure imaginary). Then, using ( |127| ), the Bogolubov coefficients can 
be read off directly from ([1321) . Using (jjgjj), W^{z) = W_ A)ff (-z), and the 
relation ^ 

|T(- + iy)\ 2 = 7r(cosh7ry) _1 



we have 

For \a ± A| >> 1, 
where 



COsh 7T 


a - X 




e 27r H cosh7r 


<r + A 


— COsh 7T 


(7- A| 



-2tt|o-+A| 



(133) 
(134) 



± A| 



(.M 2 # 2 + e 2 ^ 2 ) 1 / 2 ± e£ l + 0(# 2 O 



As mentioned above, in keeping with standard notation we have taken k > 0. 
The result for k < is obtained by changing the sign of e in the final 
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expression, since the differential equation only depends on the relative sign 
of k and e. Therefore 



exp 



2vr 



(M 2 H' 2 + e 2 E 2 ) 1 / 2 + eE ( 



(k > 0) 



\Pk 



exp 



2tt 



(M H + e 2 E 2 



2\l/2 



eEn 



(k < 0) 



(135) 



In the instanton calculation we took the particle with charge e to be the one 
to the right of the inside region, and the particle with charge — e to be the 
one to the left. From (J106|), the particle to the right has k < 0, so we have 



to use the second equation in order to compare with the instanton results. 
From (I 



and (|135|) , we have 



Therefore 



dAf 



\Pk\ 



-S E 



,dk 



-Se 



dk 



(136) 



dx " "" 2n 2n ' 

in agreement with (p. 071) . 

Apart from the agreement between these momentum distributions and 
between (|108 ) and ( p. 231 ), the time tk at which the definition of particle starts 
being meaningful [see (|129Q 1, is the same as the time of nucleation to in the 
instanton formalism [see (|106|) 1, a suggestive coincidence. 



9 Summary and conclusions 

We have computed the nucleation rates for the process of membrane creation 
by an antisymmetric tensor field in a spacetime of dimension d — N + 1, for 
N =1,2, and 3. 

To this end, we have evaluated the contribution of the relevant instan- 
tons to the semiclassical partition function. These instantons are iV-spherical 
worldsheets of radius Ro [given by fl2~T|)] embedded in a Euclidean de Sitter 
background, which is itself a <i-sphere of radius H^ 1 . The flat space in- 
stantons are obtained by taking H — > 0. We have discussed the analytic 
continuation of the instantons, describing the motion of the membranes af- 
ter nucleation. The Lorentzian solutions are spherical membranes which at 
late times expand like the scale factor in the flat inflationary FRW model. 
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The effect of de Sitter transformations on a given solution corresponds to 
space-time translations of the solution. 

To evaluate the preexponential factor of the instanton contribution, it is 
necessary to study small fluctuations of the instanton worldsheet. We have 
reviewed the covariant theory of such perturbations, according to which the 
normal displacement of the worldsheet is viewed as a scalar field <j) living 
on the unperturbed worldsheet. We have given a kinematical derivation 
of the equations of motion for 0, based on the fact that the zero modes, 
which correspond to infinitessimal translations of the instanton, have to be 
solutions. 

The evaluation of the prefactor is seen to be equivalent to the calculation 
of the effective action for a free scalar field (the field <fi mentioned above) 
living in a curved background (the iV-sphere). The functional determinants 
that arise form Gaussian integration can be explicitly calculated using (- 
function regularization. This method automatically removes all ultraviolet 
divergences. In the case N = 2, an arbitrary renormalization scale /i appears 
in the final result. Rescalings of fi are seen to be equivalent to a finite 
renormalization of the membrane tension and of a new term which has to be 
added to the action. This new term is just the Einstein-Hilbert action on the 
world-sheet, which for N = 2 is a topological invariant and hence does not 
contribute to the equations of motion. 

In flat space, the nucleation rates are obtained using the standard formula 
which relates them to the imaginary part of the free energy [see p6|)1. We 
have recovered known results for the production of kinks in 1 + 1 dimensions, 
charged pairs in 3 + 1 dimensions and 'bubble' formation in 2+1 and 3+1 
dimensions. Our results apply only to the case when the membranes are 
infinitely thin, having no internal structure. The effect of finite thickness of 



the membrane can be important in realistic field theories |2^, 25], |23], 2"?fl . We 



plan to return to this question in a forthcoming paper, where the thermal 
properties of the nucleated objects are studied |2P| . 

In de Sitter space, there is no standard procedure for the calculation of 
nucleation rates. We have introduced a heuristic prescription to obtain the 
distribution of nucleated objects during inflation, following ideas related to 
the treatment of this problem in Ref. [||]. This distribution is obtained 
from the one instanton contribution to the partition function. The param- 
eters corresponding to the zero mode rotations, which form a subgroup of 
0(d+ 1), have to be analytically continued along with the instanton so that 
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they correspond to a subgroup of the de Sitter group 0(d, 1). With these 
manipulations, the one instanton contribution to the partition function lends 
itself to interpretation as a distribution of membranes in a space of parame- 
ters. The parameters can be chosen to be the place and time of nucleation. 
In flat space this prescription reduces to the standard formula (|90f) for the 
calculation of nucleation rates. From the parameter distribution one can eas- 
ily find the size distribution of membranes in the inflationary universe, which 
turns out to be stationary and nearly scale invariant [see (P-Olj )l ■ 

For the case of pair creation, we have given the distribution in terms of 
the conserved coordinate momentum k. This distribution turns out to be 
independent of k (as expected on the grounds of scale invariance) with an 
upper cut-off k 2 hys ~ M 2 + (eE^/H 2 ), where e is the charge of the particle, 
E is the electric field, and k p h ys = e~ 2Ht k is the physical momentum. This 
cut-off corresponds to the momentum of the particle at the time of nucleation 
to (the physical momentum is subsequently redshifted). 

For the case of pairs, the calculations can be repeated using a completely 
different approach. In 1+1 dimensional de Sitter space, the Klein-Gordon 
equation for a charged scalar field coupled to a constant electric field can be 
solved in terms of Whittaker functions. We find the Bogoliubov transforma- 
tions between an 'in' state which is analogous to the Bunch-Davies vacuum 
and an 'out' state which corresponds to the definition of particles at late 
times. This definition of particles at late times is natural and unambiguous 
when the effective mass M 2 + (e 2 E 2 /H 2 ) is much larger than the expan- 
sion rate (this is essentially why in everyday experiments one does not care 
whether we live in flat space or in a de Sitter space of very small H). In this 
limit, the results can be compared with those obtained using the instanton 
formalism. Both methods agree, not only in the exponential dependence but 
also in the prefactors. This is true even when the size of the instantons is 
comparable to the horizon size. 

To summarize, our results seem to indicate that the spontaneous nucle- 
ation of defects during inflation and the decay of false or true vacuum through 
nucleation of true or false vacuum bubbles is well described by the instanton 
formalism, at least in 1+1 dimensions. Also, that our prescription for finding 
the equilibrium distribution of nucleated objects from the semiclassical par- 
tition function is correct. Although we have not attempted to give a rigorous 
justification to this prescription, we hope that with the examination of fur- 
ther examples a clearer picture will emerge. After this paper was completed, 
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we became aware of Refs. , dealing with the semiclassical approximation 
to the path integral. The methods developed in these references may be use- 
ful to provide a rigorous foundation to our prescription. Work along these 
lines is currently under way. 

Finally, some comments on negative modes. Coleman has shown ]39[] that 
in flat space and at zero temperature an instanton describing the decay of 
a metastable state has one and only one negative mode. As noted in [|], 
the instantons describing the spontaneous nucleation of defects during infla- 
tion can have more than one negative mode. Indeed, we have seen that the 
scalar fields <f>w describing the transverse displacements of the worldsheet are 
tachyonic, and that as a result each field has a negative mode with L = 0. 
Therefore, the number of negative modes for such instantons is equal to the 
co-dimension of the worldsheet. It was shown in |5j that this is not in direct 
contradiction with Coleman's theorem, which does not apply in de Sitter 
space (or in flat space at finite temperature). Still, the question remains of 
whether the wrong number of negative modes renders the instanton 'unphys- 
ical'. From the analysis of the previous Section, we think that this is not 
the case. Using the method of Bogoliubov transformations, we find that par- 
ticles are produced in de Sitter space for arbitrary co-dimensionality of the 
world-line, and that the results are always in agreement with the instanton 
results. Also, a nice feature of the prescription given in Section 7 is that the 
distribution dAf is always real regardless of the co-dimension, because the 
extra factors of % coming from additional negative modes are compensated 
by imaginary factors coming from the complexification of additional boost 
zero modes. 
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Appendix 



In this Appendix we evaluate the functional determinants for a free scalar 
field on the A^-sphere, for AT = 1, 2 and 3. As explained in section [| this 
scalar field represents small fluctuations of the worldsheet of the instanton, 
and the determinant gives the pre-exponential factor in the semiclasical par- 
tition function. The basic idea is to define a generalized zeta function (179) 
and then express the determinant in terms of £(0) and C'(0) = d(/dz\ z = 
[see Eq. (T7S1)]. For the evaluation of ((z) we follow the method of Ref.pl]. 



where the case N = A was studied (for fields of arbitrary spin). 
Al Determinant on the circle 



This is the case N — 1. From ( fH] ) and (77 ) the eigenvalues are 

A L = L 2 + M 2 R 2 = L 2 + x 2 , (L = 0, oo). 
For L ^ the degeneracy is 2 and for L = it is 1, so 



C(*) = + £ 22T* [1 + ^ 



2 \ — z 
X X 



Expanding the binomial term in powers of x/L 



Hr-f>(^r- <ai) 

we have 

oo 

C(z) = x~ 2z + £ 2c k x 2k (-l) k ( R (2z + 2k), (A2) 

fc=0 

were (r(z) = Y^l=i L~ z is the usual Riemann's zeta function. Of the coeffi- 
cients Cfc all we need to know is that 



c = 1, (A3) 



Ck = - + 0(z 2 ). (k>l) 

Then it is clear that 

C(0) = 1 + 2^(0) = 0, 
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since Cr(0) = -1/2. 

To evaluate ( ff9|) we also need C'(0)- Expanding (|A2|) in powers of z we 
have 

oo 

C(z) = -2z \n(2nx) + £ 2-x 2k ( R (2k)(-l) k + O^ 2 ), 

fc=i K 

where we have used ( R (0) = —1/2,^(0) = -(1/2) ln(27r). As a result, 

OO r> 

C'(0) = - ln(27rx) 2 + J2(-l) k T x 2k C R (2k). (A4) 

fc=i fc 

Now we need a technique to sum the series. 
For convenience one introduces the notation 

oo 

C R (z,a) = ]T L~ z , 

k=a 



so that (r(z, 1) = (r(z). Using [pC 



-l) n+1 n!G*(n + 1, a) 



one easily arrives at the expressions 



OO -| 

]T C*(2n + 1, a)z 2n = --Ma + z) + *(q - «)], (A5) 

n=0 ^ 



and 



2, 



£ Cn(2n, a)^ 2 ™ = -Ma + z) - v&( a - (A6) 

n=l ^ 

where = dlnr(z)/ciz is the digamma function. 

Now for the evaluation of ( |A4|) . Differentiating with respect to x and 
using ( |A6|) we have 

d 2 

— C'(0) = — + 2i[tf (1 + ix) - - ix)]. 
Upon integration we obtain the result 

C'(0) = -21n(2sinhvrx). (A7) 
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In the last step, the constant of integration is chosen so that ( |A7| ) agrees with 
( |A4| ) when x — > 0. Restoring x = MR we obtain the desired result fl80D . 

Note that, since £(0) = the determinant does not depend on the ar- 
bitrary parameter \x. This was to be expected, since after all for N = 1 
the path integral is just equivalent to the path integral for a non-relativistic 
quantum mechanical oscillator at finite temperature. That is, (|76|) is equiva- 
lent to (H) with V(x) = M 2 x 2 /2 and (3 = 2ttR . In (||) the right hand side 
is shorthand for 



1 (^j+l ~~ x i 



+ V{ Xi ) 



(A8) 



where e = (3/n and x n+ \ = x±. Here, there is no arbitrarity in the definition 
of the measure of integration, and so we do not expect any arbitrary scales to 
show up in the final result. The integral ( |A8| ) can be done using the methods 
of Ref. |41[] applied to the case of periodic boundary conditions, and one 
recovers the result 



A2 Determinant on the 2-sphere 

In this case the eigenvalues are given by 



A L = L(L + 1) + M Rq = (L + - + u)(L + - 



u 



u 



rr 



where u 2 = (1/4) — M 2 Rf ) . The degeneracies are given by gi = 2L + 1. It 
follows that 

oo - „,2" 

C(z) = E ^~ 2z 

n=l/2 

where n runs over the positive half integers. Expanding the binomial term 
in powers of (u/n) one has 

OO 1 

az) = Y,2c k u 2k ( R (2z + 2k-l,-). 

k=0 1 

Expanding in the neighborhood of z = we obtain 

C(*) = + u 2 + z[A(' R (-l, 1 -) + Q]+ 0(z 2 ) 



(A9) 
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where 

3 = £ ^Ch(2* - 1,^) - 2« a * (i 



Here we have used 



and the relation 



CM-l,a) = + 2 a ~~ 12' 



It is clear that £(0) = (1/12) + u 2 



To evaluate C'(0) we nee d to find Q. Using (|A5|) one easily arrives at 

After a bit of algebra 

Q = -Z7T - 3 in f _ - y \ + C + 0{2u - 3), 

where C is a numerical constant of order unity. The last term, indicated as 
0{2u — 3) vanishes when u — > 3/2, i.e., when M 2 — > — 2i?g 2 , the case we are 
interested in. 

To summarize, we have 

C(0) = i + «. 2 

('(0) = -Z7T - 3 ln(- - u) + C + 0(2w - 3), 

where we have absorbed the term 4£^(— 1, 1/2) in a new constant C". Using 
fl89|) one finds, in the limit u — > 3/2, 

(det't^o)- 2 ^])" 172 = (vR ) 7/3 Ro 3 . (A10) 

Some numerical constants have been absorbed in a redefinition of the renor- 
malization scale \x. 



52 



A3 Determinant on the 3-sphere 

In this case, from ( f44|) and ([T7]), with y 2 

K L = {L+lf-y\ 

The degeneracy is given by gi = (L + l) 2 , so 

,2 



1 - R 2 M 2 , 



C(z) = E 



n 



2-2z 



n=l 



1 - 



y 



rr 



J2cky 2k ( R (-2 + 2z + 2k), 

k=0 



(All) 



where c& are given by (|A3|) . Expanding around z = we have 
C(z) = Ch(-2) + «[2&(-2) + Q] + 0(z 2 ), 

where 



oo 2fc 

o = E V^( 2A: - 2 ) 
fc=i K 



(A12) 



It is clear that C(0) = Cfl(~2) = 0. 

Also, C'(0) = 2(' R (-2) + Q. To evaluate Q we first diferentiate (|AT2|) and 
then use ( |A6| ) to find, after some algebra, 

dQ 2 d n / • m 



Integrating, 



2 r n y 

Q = — y 2 ln(sin iry) H — -J xln(sinx)<ix 



(A13) 



where the constant of integration is chosen so that Q(y = 0) = 0, in order to 
agree with (|A12 ). 

The integral in ( |A13| ) cannot be done analytically for arbitrary ny. How- 
ever, this is not a problem, since we want to calculate the determinant only 
for M 2 — > — 3Rq 2 , which implies y = 2. Then the second term in ( |A13| ) is a 
definite integral which can be found in the tables fl37j, and we have 



Q( y ) = -y 2 ln(smny) - 4 ln2 + 3m + 0(y - 2). 
Putting it all together we have 

C(0) = 0, (A14) 

C'(0) = 2&(-2) - 4 In 2 + 3m - y 2 \n(sirmy) + 0(y - 2), 
where y = (1 - R 2 M 2 ) 1 / 2 . 
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Figure caption 



• Fig.l Euclidean de Sitter space is a rf-sphere of radius H~ l . The in- 
stantons for membrane creation and spontaneous nucleation of defects 
can be seen as spherical worldsheets of dimension N and radius Rq 
embedded in the rf-sphere. For the case of spontaneous nucleation, the 
instantons have maximal radius i? = H' 1 . For the case of membrane 
creation, the co-dimension of the worldsheet is one, d = N + 1, and the 
instanton can be obtained by intersecting the (i-sphere with a hyper- 
plane at a distance ujq from the origin. The worldsheet of the membrane 
divides the <i-sphere into two regions, which we conventionally denote 
as the inside and the outside of the membrane. The value of the elec- 
tric field in the outside region is taken to be equal to the background 
electric field before nucleation E . By Gauss' law, the electric field in 
the inside region is E — e. 
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